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Abstract 

Using spin-wave theory, we show that geometric frustration fails to preserve a two-dimensional 
spin fluid. Even though frustration can remove the inter layer coupling in the ground-state of a 
classical antiferromagnet, spin layers innevitably develop a quantum-mechanical coupling via the 
mechanism of "order from disorder". We show how the order from disorder coupling mechanism 
can be viewed as a result of magnon pair tunneling, a process closely analogous to pair tunneling in 
the Josephson effect. In the spin system, the Josephson coupling manifests itself as a biquadratic 
spin coupling between layers, and for quantum spins, these coupling terms are as large as the 
in-plane coupling. An alternative mechanism for decoupling spin layers occurs in classical XY 
models in which decoupled "sliding phases" of spin fluid can form in certain finely tuned conditions. 
Unfortunately, these finely tuned situations appear equally susceptible to the strong-coupling effects 
of quantum tunneling, forcing us to conclude that in general, geometric frustration cannot preserve 
a two-dimensional spin fluid. 

PACS numbers: 71.27.+a, 75.30.Ds 
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I. INTRODUCTION 



This study is motivated by recent theories of heavy electron systems tuned to an an- 
tiferromagnetic quantum critical point- 1,2 which propose that the formation of magneti- 
cally decoupled layers of spins plays a central role in the departure from Fermi liquid 
behavior. A wide variety of heavy electron materials develop logarithmically divergent 
specific heat coefficients and quasi-linear resistivities in the vicinity of quantum critical 
point o 1 ^! 4 ! 8 ^ 6 ^^ 18 ! 19 ^ 20 ! 21 ! 22 ^ 3 ! 24 ^ 5 ' 26 ^ 27 . Several theories explaining these unusual properties 
have been propose d 1 ! 2 ! 4 ! 5 ! 6 ! 7 ! 8 ! 18 ! 20 ' 26 . The standard model for these quantum phase transi- 
tions, proposed by Hertz and Moriya, involves a soft, antiferromagnetic mode coupled to 
a Fermi surface. Hertz-Moriya SDW theory can account for the logarithmically divergent 
specific heat coefficients and quasi-linear resistivities 1,3 , but only if the spin fluctuations 
are quasi-two-dimensional. An alternative local quantum critical description, based on the 
extended dynamical mean field theory, also requires a quasi-two-dimensional spin fluid 2 ". 
Each of these theories can only account for the anomalies of quantum critical heavy electron 
materials if the spin fluctuations of these systems are quasi-two-dimensional 1 ^ 4 ! 5 ! 6 ! 7 ! 8 . 

The hypothesis that heavy electrons involve decoupled layers of spins motivates a search 
for a mechanism that might preserve quasi-two-dimensionality in a diverse set of heavy 
fermion materials. One such frequently cited mechanism is geometric frustration^. Here, 
the idea is that frustration, naturally induced by the structure of the crystal, decouples 
layers of spins within the material^* 4, (see Fig. QJ. In this paper, using the Heisenberg 
antiferromagnet as a simple example to explore this line of reasoning, we show with the 
help of spin-wave theory that in general, zero-point fluctuations of the spin overcome the 
frustration and generate a strong interlayer coupling via the mechanism of "order from 
disorder" ^i -. 

To illustrate the main points of our argument, consider two separate layers of Heisenberg 
spins. At T = each layer is antiferromagnetically ordered, and spin waves run along the 
layers. Now consider the effect of a small frustrated interlayer coupling. In a system of 
classical spins, the layers remain decoupled in the classical ground state, and their spins 
may be rotated independently. The long-wavelength spin waves continue to run along the 
layers, and the spin fluid is quasi-two-dimensional at long wavelengths. 

In the quantum-mechanical picture, even a small interlayer coupling enables magnons 
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FIG. 1: Lattice Structure. 

to virtually tunnel between layers. An antiferromagnet can be regarded as a long-range 
RVB stated, so individual magnon transfer is energetically unfavorable, and the transfer of 
magnons between the layers tends to occur in pairs, as in Josephson tunneling (see Fig. 2). 
Interlayer magnon pair tunneling is ubiquitous in three-dimensional spin systems, frustrated 
and unfrustrated alike. So unless the interlayer coupling constant is set exactly to zero, 
magnons travel between the layers, producing a coupling closely analogous to Josephson 
coupling of superconducting layers. Such a coupling is an alternative way of viewing the 
phenomenon of "order from disorder"^i£, whereby the free energy of zero-point or thermal 
fluctuations depends on the relative orientation of the classical magnetization. 

If we use the analogy between superconductors and antiferromagnets, then spin rotations 
of an antiferromagnet map onto gauge transformations of the electron phase in a supercon- 
ductor. In a superconducting tunnel junction, the Josephson energy is determined by the 
product of the order parameters in the two layers, i.e. 

OC COS(0 2 - 4>l) 

where t± is the tunneling matrix element, A the superconducting gap energy and ipi a (/ = 
1, 2) represents an electron field in lead one and two. By analogy, in a corresponding "spin 
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junction", the coupling energy is determined by the product of the spin-pair amplitudes. 
Suppose for simplicity that the system is an easy-plane XY magnet, then 

A£ 5 (A0) ~ -^Re [(S+(i)S+(j))(S{(i')S{(f)] oc cos(2A0), 

where represents the spin raising, or lowering operator at site i in plane I, parallel to 
the local magnetization. The factor 2A0 arises because the spin-pair carries a phase which 
is twice the angular displacement of the magnetization (S + = S x + iS y ~ Se 1 ^). In other 
words, 

J2 

A£s(A0) -ry cos 2 (A0) + const, 

II 

so the interlayer coupling induced by spin tunneling is expected to be biquadratic in the 
relative angle between the spins. Clearly, this is a much oversimplified argument. We 
need to take account of the 0(3), rather than the U(l) symmetry of a Heisenberg system. 
Nevertheless, this simple argument captures the spirit of the coupling between spin layers, 
as we shall now see in a more detailed calculation. 



II. SPIN- WAVE SPECTRUM FOR DECOUPLED LAYERS 

Consider a Heisenberg model with nearest-neighbor antiferromagnetic interaction in its 
ground state defined on the body-centered tetragonal lattice. This choice of model is moti- 
vated by the structure of CePdiSii, one of the compounds for which the idea of quasi-two- 
dimensionality was originally proposed^. In this lattice structure (Fig. [TJ), square lattices 
stack with a shift of (|, |) between adjacent layers (a is the lattice constant within the 
layer). For simplicity, the distance between the layers is also a. The spins of the nearest 
neighbors in each layer are anti-parallel. In the classical ground state the spins in different 
layers are decoupled and may assume any relative alignment. 

For simplicity, let us consider just two adjacent layers, the argument being easily gener- 
alized to an infinite number of layers. The Hamiltonian is then 

H = H + V, (1) 

with 



H = H^+H^ T \ 
4 



(2) 





FIG. 2: Contrasting (a) Josephson tunneling between paired superconductors and (b) magnon 
tunneling between antiferromagnets, viewed within a resonating valence bond (RVB) picture. 

where and are the Hamiltonians for the top and bottom layers, and V is the 

interlayer coupling. 

i/„ = ^"E( s . lB,s ffi+ s Ss. ( 3 + A), (3) 

i,A 

v- = ^£(s! B, s!3 + srs£_ a ). (4) 

i,A 

Here S; B ^ (Sp^) is the spin variable defined at the site i in the bottom (top) layer. The 
vector A denotes a displacement to the nearest neighbor sites within the plane, A = (a, 0) 
or (0,a). 5 = (a/2, a/2) defines a shift between layers. Since the coupling between layers 
is small (J 1 - <C J"), we may treat this model using perturbation theory where the ratio of 
coupling constants J^/J" is taken as a small parameter. 

For our purposes, it is sufficient to consider a simple case with the spins lying in the 
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planes of the 2-dimensional lattice. At sites i = (la, ma) and i + 8 — (la+ |a, ma + ^a) the 
spins are 



= S(-l) l+m , SP B) = 0; (5) 
Sf + ( P = S(-iy +m+1 cos0, S^P = 5(-l)' +m+1 sin0; (6) 

where X and Y are mutually perpendicular directions in the plane, I and m are integers. 

Following a standard procedur o 12 i 13 , we use the Holstein-Primakoff approximation for 
the spin operators to determine the spin-wave spectrum. The single-layer Hamiltonian 
becomes 

#(«) = —ANS 2 J^ + [sSJ^a+^a*^ + SJ 11 (q) [o^o£J + /i.e.]] , (7) 

q 

and on diagonalization the Hamiltonian Hq for the decoupled layers can be written as 



a=T,B q 

The ground state energy of the decoupled two-layer system is then 



E = -8NS(S +1)J" + (9) 

q 

u\ defines the spectrum of spin waves propagating in each of the layers 



u><\ = AS J" J 4 — [cos q^a + cos q^a} 2 . (10) 



III. MAGNON PAIR TUNNELING BETWEEN THE LAYERS 



Now we express the perturbation V in (0J) in terms of 



+(df) (a) 
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where and are defined as 



K = 2J± 



qt.a q v a . q^a . q v a 
cos(-^— J cos(-^-J + sm(-^-J sm(-^-J cos < 

Z Zi Z Z 



(12) 



Bi = 2J 1 - 



cos(-|-) cos(^) - sm(-|-) sm(-^-) cos0 



In terms of bq^ a \ bq^ 



(13) 



V = SJ2 + fc+^O + V ph , (14) 

q 

where a q = (A^ sinh 2w q + B^ cosh 2« q ) describes the amplitude for magnon pair tunneling, 



and V ph = S'^ q (A q - cosh2w q + B^ sinh 2w q )(& q (T ' ) & q B ' ) + bq (n> bq ') describes single magnon 
tunneling between layers. 

It is straightforward to see that the particle-hole terms do not affect the ground-state 
energy, for V ph \GS >= 0, where \GS > denotes the ground state wave-function of the 
system. The second order correction to the ground state energy E is then 



AE, 



(2) 



E 



< \\V\GS > 



E\ — E, 



GS 



(15) 



where |A > denotes a state with two magnons being transfered between layers. Thus, 



AE^ = -J2S 2 al/(2ul). 



(16) 



To understand the nature of coupling between the layers (dipolar or quadrupolar) , let us 

(2) 

retrive the dependence of AEq on the angle 0. 



AE. 



(2) 



S (J r ) 2 
'2 JW 



[C + C 2 cos 2 0] 



2 jr [(Co + % + ^c s24 



(17) 



The particular form of the coefficients is 
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T , 

Co = J ^^|y 1 ~ ^[cosx + cosy] 2 cos 2 I cos 2 |(l - ^[cosx + cosy]) 2 , (18) 



^ 2 = J (2 ) 2 \f ~ 4 [ cosx cos ^ 2 cos2 ^ s ^ n2 2 sm2 2 2 ^ CQS X CQS ^]) 2 - (1^) 

— 7T 

The interlayer coupling is indeed quadrupolar in nature, as foreseen earlier. Moreover, 
there is no small parameter, and for small S, when J 1 - ~ j", this coupling is not weak. 



IV. DISCUSSION 



In the above calculation, we considered an ordered Heisenberg antiferromagnet at zero 
temperature. In practice, provided the spin-spin correlation length £ is large compared with 
the lattice constant a, £ ^> a, a biquadratic interlayer coupling will still develop. More- 
over, at finite temperatures, thermal fluctuations will produce further interlayer coupling. 
Both thermal and quantum interlayer coupling processes are manifestations of "order from 
disorder". The main difference between the thermal and quantum coupling processes lies 
in the replacement of the magnon occupation numbers with a Bose-Einstein distribution 
function, and in general both the sign and the angular dependences of the two couplings 
are expected to be the same 9 . In general, geometrical frustration is an extremely fragile 
mechanism for decoupling spin layers and will always be overcome by quantum and thermal 
fluctuations. Our work was motivated by heavy electron systems. These are much more 
complex systems than insulating antiferromagnets, but if our mechanism for the formation 
of two-dimensional spin fluid is to be frustration, it is difficult to see how similar interlayer 
coupling effects might be avoided. We are led to conclude that for the hypothesis of the 
reduced dimensionality of the spin fluid in heavy fermion materials to hold, a completely 
different decoupling mechanism must be at work. 

In the special case of XY magnetism there is, in fact, one such alternative mechanism, 
related to "sliding phases". Some heavy fermion systems, such as YbRh 2 Si2, are XY-like, 
most others, such as CeCu^, are Ising-like. It is, therefore, instructive to consider whether 
the sliding phase mechanism might be generalized to Heisenberg or Ising spin systems to 
provide an escape from the fluctuation coupling that we have discussed. 
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The existence of a " sliding phase" in weakly coupled stacks of two-dimensional (2D) XY 
models was predicted by O'Hern, Lubensky and Toner—. In addition to Josephson interlayer 
couplings, these authors included higher-order gradient couplings between the layers. In the 
absence of Josephson couplings, these gradient couplings preserve the decoupled nature 
of the spin layers, only modifying the power-law exponents of the 2D correlation functions, 
(SiSj) ~ r~ v . As the temperature is raised, Josephson interlayer couplings become irrelevant 
above a particular "decoupling temperature" T d . One can always select interlayer gradient 
couplings to satisfy T d < T KT and produce a stable sliding phase in the temperature window 
T d <T< T KT . 

To see this in a little more detail, consider the continuous version of the Hamiltonian of 
two layers of XY models, H = H + V, where H is a sum of independent layer Hamiltonians 
and V is the usual Josephson-type interlayer coupling 

#o = y J rfV[V^T(r)] 2 + y J d 2 r[V ±< Mr)] 2 , (20) 

V = J L J d 2 r cos[0 T (r) - B (r)]. (21) 

At low temperature, when the interlayer coupling J 1 - is zero, the average of the intralayer 
spin-spin correlation function with respect to Hq is 

(0 2 (r)) o = r / log(L/6), (22) 

and 

(cos[0(r) - «/>(0)])o ~ (W, (23) 

where i] = T/27rj", L is the sample width and b is a short- distance cutoff in the XY plane. 

The average of Josephson interlayer coupling V scales as (V)o ~ L 2_r? , so Josephson 
couplings become irrelevant at Td = 47rj". At temperatures above the Kosterlitz-Thouless 
transition temperature T^t — vrj"/2, thermally excited vortices destroy the quasi- long- 
range order and drive the system to disorder. In this simple example, it happens that 
Td > Txt, which does not permit a sliding phase. However, higher-order gradient interlayer 
couplings between the layers, when added to this model, suppress T d below Tkt, producing 
a stable sliding phase for T d < T < Tkt- 
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So can the sliding phase concept be generalized to Heisenberg spin systems? A sliding 
phase develops in the XY model because power-law spin correlations introduce an anoma- 
lous scaling dimension, but unfortunately, a finite temperature Heisenberg model has no 
phases with power-law correlations 15 . In general, biquadratic interlayer couplings will al- 
ways remain relevant in Heisenberg models. In the quantum-mechanical picture, as soon 
as a frustrated interlayer coupling is introduced, the order-from-disorder phenomenon*^ 
generates a coupling A ~ SJ 1 - 2 / J" between the layers: 

H = i \ E( V ^) 2 + \ Eft - ^) 2 > ( 24 ) 

J i i 

where p = S 2 a 2 jK This coupling gives us a length scale Zo determined from (/o)~ 2 ~ A/p 
or l ~ a VSJ^/J 1 . Once the spin correlation length £ ~ a exp(2n S 2 /T) within a 
layer grows to become larger than l , i.e. Iq < a exp(27r J^S 2 /T), a 3D-ordering phase 
transition occurs. An estimate of the 3D-ordering transition temperature is then T c ~ 
2ttJ^S 2 / \n(y/SJ\\ I J -1 ). The answer is essentially identical in the classical picture, for here, 
thermal fluctuations generate an entropic interlayer coupling A ~ m&x(SJ ±2 / J", TS 2 ), 
so at high enough temperatures, for large S, A ~ TJ ±2 /J^ 2 , l ~ SaJ^ 3 / 2 / J ± ^/T. A 
classical estimate of the 3D-ordering temperature is T c ~ 2-kJ^S 2 / ln( / J L ). 

Another interesting question is whether XY models permit sliding phases at T = 0. The 
decoupling temperature, as found by O'Hern, Lubensky and Tonei 1 ^, is 

UP) = TZ^TZP (25) 
Jo J p 

One sees no obvious mechanism of suppressing Td to zero. A 2D sliding phase is equivalent 
to a 3D finite temperature sliding phase, so the existence of a sliding phase in the XY model 
at zero temperature would mean a power-law phase in 3D XY model. Since no power-law 
phase exists in 3D XY-like systems, sliding phases at T = are extremely unlikely. In 
conclusion, the sliding phase scenario also fails to provide a valid general mechanism for 
decoupling layers in Ising-like and Heisenberg-like systems. 

Let us return momentarily to consider the implications of these conclusions for the more 
complex case of heavy electron materials. It is clear from our discussion that simple models 
of frustration do not provide a viable mechanism for decoupling spin layers. One of the 
obvious distinctions between an insulating and a metallic antiferromagnet is the presence 
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of dissipation which acts on the spin fluctuations. The interlayer coupling we considered 
here relies on short-wavelength spin fluctuations, and these are the ones that are most 
heavily damped in a metal. Our exclusion of such effects does hold open a small possibility 
that order- from-disorder effects might be substantially weaker in a metallic antiferromagnet. 
However, if we are to take this route, then we can certainly no longer appeal to the analogy 
of the insulating antiferromagnet while discussing a possible mechanism for decoupling spin 
layers. 
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